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Abstract 

In this work we investigate small clusters of bosons using the hyper spherical harmonic basis. We 
consider systems with A = 2,3,4,5,6 particles interacting through a soft inter-particle potential. 
In order to make contact with a real system, we use an attractive gaussian potential that reproduces 
the values of the dimer binding energy and the atom-atom scattering length obtained with one of 
the most widely used 4 He- 4 He interactions, the LM2M2 potential. The intensity of the potential 
is varied in order to explore the clusters' spectra in different regions with large positive and large 
negative values of the two-body scattering length. In addition, we include a repulsive three-body 
force to reproduce the trimer binding energy. With this model, consisting in the sum of a two- and 
three-body potential, we have calculated the spectrum of the four, five and six particle systems. 
In all the region explored, we have found that these systems present two states, one deep and 
one shallow close to the A — 1 threshold. Some universal relations between the energy levels are 
extracted; in particular, we have estimated the universal ratios between thresholds of the three-, 
four-, and five-particle continuum using the two-body gaussian potential. They agree with recent 
measurements and theoretical predictions. 
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I. INTRODUCTION 



Systems composed by few atoms having large value of the two-body scattering length, a, 
with respect to the natural length, £, fixed by the atomic potential, have been the object 
of an intense investigation both from a theoretical and experimental point of view (for 
recent reviews see Refs. [l 3 ] ) . In fact, they present universal properties: for example, the 
three-body system displays the Efimov effect 0, Q], that means the appearance, in the limit 
a/ £ — > oo, of an infinite set of bound states accumulating toward the three-particle threshold; 
moreover, the three-body spectrum has a discrete-scale symmetry, with an universal ratio 
between the n-th and n + 1-th levels E^ +1 / E^ = e _2?r / so . The scaling factor depends 
only on the ratio between particle masses, and for identical bosons of mass m it reads 
e -27r/s ~ 1/515.03 (with s ps 1.00624). The finite value of £ implies the existence of a 
three-body ground state £3 whose value reflects the short range physics, and that, together 
with the discrete-scale symmetry, completely determines the spectrum. In realistic cases the 
ratio ajl is large but finite; thus, the three-body spectrum reduces to a finite number of 
states. 

A remarkable property in a — > 00 limit appears in the four-body system: two states 
E2'°, E2' 1 are attached to each trimer state EV, one deep and one shallow having universal 



-0 



ratios, E n /jEl « 4.6 and E%' /Eg « 1.001 [6|-|8(; the two lowest four-body states, E\ 
E°'° and E\ = E / , are real bound states. These properties have been studied for large 
positive and large negative values of the scattering length in the (a -1 , ft) plane, with k = 
sign(£')[|£'|/(/i 2 /m)] 1 / 2 , constructing what is normally called an Efimov plot 9|. 

There are very few studies of the spectrum of small bosonic clusters beyond A = 4. In 
addition to the specific problems related to the solution of the Schrodinger equation for 
more than four particles, the atom-atom realistic potentials present a strong repulsion at 
short distances which makes the numerical problem more difficult. Specific algorithms have 
been developed so far to solve this problem: the Faddeev equation has been opportunely 



modified 



101 ] . the Hyperspherical methods resorted either to the hyperspherical adiabatic 



(HA) expansion (for a review see Ref. [U[), or to the correlated hyperspherical harmonic 
expansion (CHH) [l^j]. However, due to the difficulties in treating the strong repulsion, few 
calculations exist for systems with more than three atom, Fo r example, in Ref. Q the 
diffusion Monte Carlo method has been used to describe the ground state of 4 He molecules 



up to 10 atoms, and in Ref. [l^] a very extended calculation has been done in the four 
helium atom system. On the other hand, descriptions of few-bosons systems using soft-core 
potentials are currently operated (see for example Refs. 

aa. 

The equivalence between hard- or soft-core-potential descriptions has been discussed in 
Refs. [if], [l?]], in which an attractive soft 4 He- 4 He gaussian potential has been used to 
investigate the three-atom system. The soft-two-body potential was designed to reproduce 
the helium dimer binding energy E 2 , the 4 He- 4 He scattering length a, and the effective range 



r of the LM2M2 potential 



one of the most used 4 He- 4 He interactions. In this context 



the soft gaussian potential can be considered as a regularized-two-body contact term in an 



Effective Field Theory (EFT) approximation of the LM2M2 19j; this is possible because 



of the scale separation between the 4 He- 4 He scattering length, a = 189.41 a.u., and the 
natural length I = 10.2 a.u., which is the van der Waals length calculated for the LM2M2 
potential [1]. 

In the two-body sector and in the low-energy limit, the two potentials predict similar 
phase shifts, therefore, even if their shape is completely different, they describe in an equiv- 



alent way the physical processes in that limit [19|. The equivalence is lost as the energy is 
increased, when the details of the potential become more and more important. When the 
soft interaction is used in the three-body sector, a new three-body-contact term is required 
to reproduce the ground-state-binding energy of the helium trimer given by the LM2M2 
potential. This term is introduced by means of a gaussian-hypercentral three-body force, 
whose strength is tuned to reproduce the LM2M2 ground state binding energy of the three- 



atom system. In Ref. [16| the quality of this description has been studied for different 
ranges of the three-body force by comparing the binding energy of the excited Efimov state 
and the low-energy helium-dimer phase shifts to those obtained with the LM2M2 potential. 
In Ref. |l3] the spectrum of small clusters of helium atoms has been investigated up to six 
particles maintaining however fixed the values of a and E 2 as given by the LM2M2 potential. 

In the present work we extend the analysis of the A = 3 — 6 bosonic spectrum to the 
(a~ l ,n) plane. We have modified the strength of the LM2M2 potential in order to cover 
the region of negative values of a up to aP_, with this value indicating the threshold of 
having a three-body system bound. We have also increased the intensity of the interaction 
in order to extend the analysis to positive values of a in which the universal character of 
the system starts to be questionable, i.e, when the ground-state E® approaches the natural 



energy Ei = —h 2 /ml 2 , which delimits the Efimov window. 

Associate with the different values of a of the modified LM2M2 potential, we have con- 
structed a set of attractive gaussian potentials with the strength fixed to reproduce the 
low-energy data of LM2M2. Moreover, the modifications of the LM2M2 produce different 
values of the A = 3 ground state energy E®; accordingly, we introduce a soft three-body force 
devised to reproduce those values along the (a -1 , k) plane. Within this model, consisting in 
the sum of a two- and three-body potentials, we have calculated the spectrum of the four, 
five and six particle systems. 

Two different calculations have been performed in the present work. From one side we 
have calculated the A = 3 ground state and excited state, E® and E\, using the LM2M2 
potential and its modification, in order to construct the corresponding Efimov plot. Since 
this potential present a strong short-range repulsion we have used the CHH expansion as dis- 
cussed in Ref. [l2[. One the other side, when using the soft-core potential model in systems 



with A > 3, the numerical calculations were performed by means of the non-symmetrized 



hyperspherical harmonic (NSE 



by the authors in Refs. 
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expansion method with the technique recently developed 
-22|. In this approach, the authors have used the Hyperspher- 
ical Harmonic (HH) basis, without a previous symmetrization procedure, to describe bound 
states in systems up to six particles. The method is based on a particular representation of 
the Hamiltonian matrix, as a sum of products of sparse matrices, well suited for a numerical 
implementation. Converged results for different eigenvalues, with the corresponding eigen- 
vectors belonging to different symmetries, have been obtained 22(. In the present work, 
since we are dealing with bosons, we only consider the symmetric part of the spectrum. In- 
terestingly, we have observed that in all the region explored the A = 4, 5, 6 systems present 
two states, one deep and one shallow close to the E\_ x threshold. To gain insight on the 
shallow state, for a selected value of a, we have varied the range of the three-body force 
and we have studied the effect of that variation in the A = 4,5,6 spectrum. In the range 
considered, the variation produces small changes in the eigenvalues, but they are crucial to 
determine if the shallow state is bound or not with respect to the A — 1 threshold. This 
analysis confirms, at least in one zone of the Efimov plot, previous observations that each 
Efimov state in the A = 3 system produces two bound states in the A = 4 system, and 
extends this observation to the A = 5, 6 systems. 

Finally, we have extended the calculations of the A = 4 and A = 5 systems up to the 



four- and five-particle thresholds using the simple two-body-gaussian potential; the ratios 
between the thresholds are in agreement with previous theoretical results {?], [2^ and with 



experiments 



24 



27|. 



The paper is organized as follows. In Section II we describe the two- and three-body 
forces we used in our calculations to reproduce the LM2M2 values. In Section III we discuss 
the Efimov plot for three particles. In Section IV the results for the bound states of the 
A = 3, 4, 5, 6 clusters are discussed whereas the conclusions are given in the last section. 

II. SOFT-CORE TWO- AND THREE-BODY POTENTIALS 

As mentioned in the Introduction, we use the LM2M2 4 He- 4 He potential as the reference 
interaction, with the mass parameter fixed to h 2 /m = 43.281307 (a.u.) 2 K. In order to 
explore the Efimov- (a -1 , k) plane, we have modified the LM2M2 interaction as following 

Vx(r) = X-V LM 2M2(r) . (1) 



Examples of this strategy exist in the literature [12j, |28(. We have varied A from A = 0.883, 
where a = a°_ = —43.84 a.u., up to A = 1.1 corresponding to a = 44.79 a.u., as shown in 
FigHJ The unitary limit is produced for A ~ 0.9743. When A = 1 the values of the LM2M2 
are recovered: a = 189.41 a.u., £2=-1.303 mK and r = 13.845 a.u.. 



Following Refs. Ill, ll6|, [l7j] we have constructed an attractive two-body gaussian (TBG) 
potential 

V(r) = V e~ r2 / R o , (2) 

with range Rq = 10 a.u., and we have varied the strength Vq in order to reproduce the values 
of a given by V\(r), as shown in Fig. [2j For A = 1 with the strength Vq = —1.2343566 K we 
reproduce the LM2M2 low-energy data, E 2 = —1.303 mK, a = 189.42 a.u., and r = 13.80 
a.u.. The use of the TBG potential in the three-atom system produces a ground state 
binding energy appreciable deeper than the one calculated with V\(r). For example, at 
A = 1 the LM2M2 helium trimer ground state binding energy is 126.4 mK whereas the one 
obtained using the two-body-soft-core potential in Eq. (T2J) is 151.32 mK. A smaller difference, 
though still appreciable, can be observed in the first excited state. 

In order to have a closer description to the A = 3 system obtained with the modified 
LM2M2 potential, we introduce the following (repulsive) hypercentral-three-body (H3B) 
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FIG. 1. The scattering length a, in units of £, as a function of the parameter A, calculated with 
the modified LM2M2 potential V\(r). The range of variation of A is between A = 0.883, which 
corresponds to the disappearance in the continuum of the excited three-body state, and A = 1.1. 
The unitary limit is obtained for A ~ 0.9743. 

interaction 

W(p 123 ) = Wo e~ p2 ^ /p o , (3) 

with the strength Wo tuned to reproduce the trimer energy E® obtained using V\(r) for all 
the explored values of A, as shown in Fig. |3j Here pl 23 = |(^i2 + r 23 + r 3i) i s the hyperradius 
of three particles and po gives the range of the three-body force, or, in the spirit of EFT, the 
cut-off of the three-body-contact interaction; therefore, it is not independent of -Ro ; which 
is the cut-off of the two-body-contact force, and in fact it should be po = Rq, as shown in 
and so we fixed po = 10.0 a.u. A different criterion to fix the three-body force was 
given in Ref. 23[ in which the condition p ^ r o nas been used. In this case the (repulsive) 
three-body force is used to push the trimer spectrum high in energy in order to verify as 
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FIG. 2. The strength Vq of the gaussian two-body-potential as a function of the parameter A. 
The values are tuned to reproduce the scattering length a given by the modified LM2M2 potential 
Vx(r). 

close as possible the universal ratios E^ +1 /E^ = e~ 27T ^ s ° already at n = 0. With the LM2M2 
interaction this relation is only approximate verified; for A = 1 we have E®/E\ ~ 56, whereas 
at the unitary limit E®/E\ ~ 525, very close to the universal ratio. 



III. THREE-BODY EFIMOV PLOT 

The calculations for A = 3 have been performed using the CHH expansion. Since V\(r) 
is obtained multiplying the LM2M2 potential by a global factor A, it inherits the strong 
short range repulsion; in this direct use of the HH basis to compute the bound 

states is not feasible since it would be necessary to include an enormous number of basis 



elements in the expansion 



29|. The use of the CHH expansion circumvents this problem 




by the introduction of a correlation factor of the Jastrow type. The method is described 



in Ref. 12j and it allows to achieve similar accuracy as other techniques. As an example, 
in Table [I] we show the results for the ground state E® and the excited state £3 at A = 1 
(in this case the results of the LM2M2 potential are recovered), and at the unitary limit 
(A = 0.9743). These results have been obtained using the CHH basis up to a value of the 
grand-angular momentum K = 160. 

As a byproduct of the tuning procedure of the three-body strength Wq, we have con- 



structed, as was previously done for instance in Refs. [28|, [30J, the Efimov plot shown in 
Fig. HI In the figure we report calculations of E® and E\ as functions of a done both with 
the V\{r) and the TBG potential. When the TBG+H3B potential is used, the results coin- 
cide with those of V\{r) and are not reported in the figure. In addition, we draw the dimer 
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potential 


£° (mK) E 


I (mK) 




A = 1 




V\(r) 


-126.4 


-2.27 


TBG 


-151.3 


-2.48 


TBG+H3B 


-126.4 


-2.31 


A 


= 0.9743 




V x (r) 


-83.99 


-0.16 


TBG 


-103.4 


-0.20 


TBG+H3B 


-83.99 


-0.16 



TABLE I. The ground state E^ and the excited state .Eg of the three-boson system calculated with 
the modified LM2M2 potential V\(r), the TBG potential, and the TBG potential plus the H3B 
potential at A = 1, that corresponds to the original LM2M2 potential, and in the unitary limit, 
A = 0.9743. 



energy E 2 , calculated using the V\(r) potential. In order to show these quantities together 
in the figure we have used the fourth root of the energy (in units of Eg) as a function of the 
square root of a -1 (in units of £). In the region analyzed, the results are inside the Efimov 
window; in fact, the scattering length is still much larger than the natural length £, and the 
ground-state energy £3 is above the natural value Eg. 

Looking in Fig.HJat negative values of a, it is possible to identify the value of the scattering 
length a}_ at which the excited state E\ disappears. For V\(r), this value is a}_ ~ — 975 a. u., 
whereas using the TBG potential it results a\ ~ — 752 a. u.. The next interesting point 
appears at a°_, when the three-body cluster is no more bound, so that E® approaches zero. 
Using V\(r) this happens at oP_ ~ —48.1 a.u., whereas using the TBG potential alone it is 
a°_ « -43.3 a.u.. 

The ratio a Q _/a l _ has been predicted to have an universal value (a°_ /a_)thcory = ^^•^ fl' m 
the only experiment which measures the two thresholds, Ref. 241, the ratio is (a la 1 ) 

J L 1 * — ' 1 \ —1 —/ experiment 



21.1. In our case we obtain (a°_/a\.) TBG = 17.4, and (a° / a -)v x M = 20.3, which is closer 
to both theoretical and experimental values. The absolute position of a°_ is not predicted 
by the theory of Efimov physics and, in that sense it can be considered as not an universal 
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FIG. 4. (color online). Efimov plot for A = 3. We report the ground and excited state energies, E® 
and -E3, in units of Eg, as & function of £/a, both for the modified LM2M2 potential V\(r), and for 
the TBG potential. Following the literature, we really draw the fourth root of the scaled energies 
as a function of the square root of the scaled-inverse scattering length; using this trick, the ratio 
between excited and ground energies is greatly reduced, allowing for the graphical representation 
of both curves on the same scale. We also report the A = 2 binding energy. 

quantity; however, it has been the subject of experimental measurements which give more 
or less the same value in units of mean scattering length a = 0.955978 i/2 for different 



atoms, a°_ = —(9 ± 1) a 3l| • In the present calculations we obtain (g^_)tbg — —8.9 a and 



(a°-.)v x (r) = -9.9 a. 

In addition, we discuss the universal character of the shallow state E\. Using the Efimov's 
radial law j^j it is possible to obtain an equation for this trimer binding energy as a function 
of a. It reads 

El + — 2 =exp[A(0/s }^, (4) 
ma 2 m 
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where is the wave number corresponding to the energy El = h 2 n1/m = 0.156 mK 
at the resonant limit and tan£ = —{mE\/h 2 ) l l 2 a. The function A(£) is universal and a 
parametrization in the range [— n, — 7r/4] is given in Ref. [l|. It verifies A(— vr/2) = and, 
from the very precise result A(— 7r/4) = 6.02730678199 and A(— tt) ~ —0.89 Q, it is possible 
to determine the values a* (at which El = E 2 ) and a}_ (at which El — 0). In order to analyze 
the universal character of the calculated energies El using the TBG and TBG+H3B potential 
models, in Fig. [5] we compare them to the values of Eq.(j4]). By construction the energies El 
at the resonant limit coincide with those of Eq.(j4]). It is possible to see that the calculated 
energies using the TBG potential, Fig.[5]upper panel, and TBG+H3B potential, Fig. [Slower 
panel, reproduce the universal behaviour close to the resonant limit. The small differences 
observed at finite values of a, especially close to the critical values a* and a^_, are due to 
effective-range corrections which are automatically included in our approach. Moreover, E\ 
does not disappear in the atom-dimer continuum at a* but follows very close the E 2 curve 
from below. 

To conclude, we further analyze the universality looking at the correlations between the 
three-body ground and excited states, as has been proposed in Ref. 32J. In Fig. [6] we trace 
the square root of the excited-trimer energy, measured from the two-body dimer, in units of 
the trimer-ground state energy, as a function of the dimer energy, always in units of trimer- 
ground state energy. The Efimov's universal-radial law Eq. (j4j) gives the universal curve in 
this plot; we see that as far as the dimer is very shallow, the calculated points are very close 
to the universal curve. They depart from it when corrections due both to finite scattering 
length and to non-zero effective range become sizeable. This non-universal effect is more 
important for the TBG case, probably due to the lack of the three-body corrections. 



IV. EFIMOV PLOT FOR A = 4, 5, 6 CLUSTERS 

The calculations for the A > 3 systems are performed using NSHH basis. The method 
been recently used to describe up to six nucleons interacting through a central potential 



nas 
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331 ] and six bosons using a two-body plus a three-body force 17j. The Hamiltonian 



matrix is obtained using the following orthonormal basis 



( P n\m[K)) = {^ +1)/ \lj^^Lt\M^ p/2 )y^ (5) 
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FIG. 5. (color online). Comparison between the excited three-body energy E% and the theoretical 
universal value given by Eq. (j3J), using the TBG potential (upper panel), and the modified LM2M2 
potential V\(r) (lower panel). The two-body eneggy E<i calculated with V\(r) is also shown. The 
calculated and theoretical curves agree around the resonant limit, and the differences close to 
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FIG. 6. (color online). Correlations between the ground and excited £"3 states of the trimer. 
We compare the universal correlation, obtained by means of Eq. (jlj), to the calculations made 
using the full V\(r) potential and the TBG potential. The agreement is good close to the unitary 
point, where the dimer energy E2 is small. The deviations become significant when the finite 
effective-range effects become non-negligible. 

where Lm\(3p) is a Laguerre polynomial with a = 3N — 1 (N = A — 1) and (3 a variational 
non-linear parameter. The function y^f{VL^) are the HH functions with grand-angular mo- 
mentum K, and total angular momenta L and magnetic number M. The Hamiltonian matrix 
is not constructed, but using properties of HH is expressed as an algebraic combination of 
sparse matrices, allowing for an efficient research of the lowest eigenvectors/eigenvalues. A 
full discussion of the NSHH method is given in Refs. 17|, |22 |. 

After solving the A = 3 problem for bound states, used to fix the strength of the H3B 
force, we have diagonalized the Hamiltonian for A = 4,5,6 bodies using the TBG and 
TBG+H3B potentials. The results are given in Fig. [7] in two scaled-(a _1 , k) plots, one 
obtained with the two-body potential alone (upper panel) and one with the two-body plus 
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three-body interactions (lower panel). In the first case, with only the TBG potential, we 
observe that the spectrum of the systems A = 4, 5, 6 presents two bound states, one deep 
and one shallow, for all values of a studied. When the repulsive three-body force is included, 
the spectrum moves up and we can observe that the excited state E A disappears for A = 5, 6 
for negative values of the scattering length as a approaches a l _. This fact is better shown in 
Fig. [8j where the differences E A — E l A+1 have been plotted as a functions of Ija. Whereas 
the differences E\ — E\ and E® — El are positive along the whole range, indicating that 
the states E\ and E\ are bound, the differences E\ — E\ and E\ — E\ result negatives as a 
goes to the negative region, so at some value of a the excited states E\ and E\ are no more 
bound. The determination of the point where the transition happens can be determined by 
looking at the convergence of the states E\ and E\, as can be seen in Table HT1 where we 
report the convergence pattern using the TBG+H3B potential for A = 4, 5, 6 at the point 
A = 0.9 where both E\ and E\ are not bound. They remain above the E® and E® threshold 
respectively. In the table we have shown the three maximum values of K considered in the 
present calculations. 

Moreover, the fact that the states E\ and E\ are bound or not depends also on the range 
of the three-body force po- I n order to analyze this relation, we have varied po aX \ = 0.9 as 
well as at the unitary limit. For each value of po the strength of the three-body potential has 
been fixed to reproduce the trimer binding energy El as before. The results for A = 4, 5, 6 
at A = 0.9 are shown in Fig. [9j As can be seen, the excited states are recovered as bound 
states for values of p ~ 18a.u.. 

To make contact with the analysis of Ref. 34] we have calculated \/\E\ — E® \ / \ E% \ = 



0.070 and a/ | E® | / 1 E\ \ = 0.434 at the unitary limit. These two values correspond to a point 
in the plot given in Fig.l of that reference lying very close to line giving the relation of these 
two quantities at the unitary limit. In addition, in Fig. [TD]we analyze the relation between 
E\ +1 (upper panel) and E\ +l (lower panel) with E A and E\ respectively, as a function of 
the scattering length for A = 4, 5, 6. We can observe a linear dependence in all cases except 
for a small curvature in the E® vs. E® and E\ vs. E\ curves close to the point in which E® 
goes to zero. These curves display the universal character of these clusters as their spectrum 
is determined by two parameters, a and E®. 

Besides, the universal ratios E A /E® and E A+l /E A can be studied at A = 1. They are: 
EH 'El = 4.5, EH 'El = 10.4, El/El = 18.4 and E\/ El = 1.020, E\/ E Q 4 = 1.009, E%/E% = 
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1.016. These ratios are in close agreement to those obtained in the literature 23J. At the 
unitary limit the ratios E A /E® move a little bit from the universal values showing some 
dependence on the form of the soft potential whereas the ratios E A+1 /E A show stability. 
At A = 0.9743 they are: E%/E$ = 5.3, E%/E$ = 13.0, E%/E$ = 23 .4 and E\j ' E\ = 1.026, 
El/ El = 1.004, E\/E% = 1.006. 

Finally, using the TBG potential, we have extended the calculations for A = 4 and A = 5 
systems up to the four- and five-particle thresholds in order to calculate the ratios between 
the different thresholds and to compare our results to previous calculations and experimental 
outcomes. Our results are summarized in Fig. [TTJ Denoting with aj° and a_! the four- 



particle thresholds of the ground and excited state respectively, we have a 



1,1 



-39.8 a.u., 



and aj° ~ 



19.6 a.u.; equivalently, with respect to the three-particle threshold a°_ we get 
(a_! / a°)TBG ~ 0.92 for the excited state, and (a 4 _l°/a°_ )tbg ~ 0.45 for the ground states. 



These values agree very wel 

□ 

predicted by the theory p, 



with what measured in the experiments 24j-|26| and with what 



35] The five-particle thresholds read a 



5.0 



■12.5 a.u. for 



the ground state and a!! 1 ~ —18.7 a.u. for the excited state. Equivalently, the ratios with 
respect the four-particle threshold are (a J /qJ°)t bg ~ 0.64 and (o^/^-^tbg »s 0.95, still 



in agreement with previous theoretical prediction 



23( and with recent experiments 27]. 



V. CONCLUSIONS 

In the present paper we have discussed the spectrum of bosonic systems up to six particles 
interacting through a two-body potential having a large two-body scattering length (with 
respect to the effective range). The three-body scale has been fixed using a scaled Helium- 
Helium potential. The scope was to extend previous studies on the Efimov physics done 
in the three- and four-body systems. We have observed that, similarly to the four-body 
system, the five- and six-body systems present two bound states, one deep and one shallow. 
It seems that this type of spectrum is to some extend universal depending only on the 
condition a ^> tq in the two-body system. This condition produces a geometrical series of 
bound states in the three-body system and, attached to each of these states, a two-level 



spectrum has been showed to appear in the four-body system js _ S|- However they are true 
bound states only in correspondence to the lowest trimer level. The other states appear 
as resonances embedded in the continuum of four particles. It is possible to introduce a 
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FIG. 7. (color online). Energies of the A = 3 — 6 ground and excited states, E^, E\, as a function 
of a" 1 , using the two-body gaussian potential (upper panel), and using the two-body plus the 
hypercentral three-body force (lower panel). In both panels we also give the two-body ground- 
state energy E^ calculated with the LM2M2 potential. 
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FIG. 8. The difference AE = E\ - E\ +1 for the indicated cases as a function of the inverse of a 
for the TBG+H3B potential. The particular cases at A = 1 and 0.9 as well as at the unitary limit 
are indicated as vertical lines. 

repulsive three-body force that eliminate all the trimer states below one specific level. In 
this way, the two level spectrum of the four-body system attached to this trimer ground state 
will become true bound states. Also the universal character of the spectrum will be more 
evident as the repulsive three-body force will push more and more the particles faraway. 

Though the analysis of the bosonic spectrum can follow the strategy illustrated above, in 
the present paper we follow a different one based on the physics introduced by the two-body 
potential. In nuclear systems as well as in many atomic systems the two-body interaction 
has a sharp repulsion at short range followed by a very weak attractive part that produces 
very shallow dimers as for example the deuteron or the two-helium molecule. The particles 
are located in the asymptotic region and do not feel the details of the interaction. There- 
fore, we can introduce a soft potential to be considered as a regularized-two-body contact 
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Po (a.u.) 



FIG. 9. The difference AE = E° A -E\ +1 at A = 0.9 as a function of p for the TBG+H3B potential. 



term in an EFT approximation of the original potential [19|. In the three-body sector, a 
three-body-contact term is required to reproduce the ground-state-binding energy of three- 
particles introduced here by means of a gaussian-hypercentral three-body force. Using this 
potential model we have calculated the three-body spectrum and we have analyzed its uni- 
versal character comparing the energy of the shallow state E\ to the Efimov's equation 
for the binding energy. Furthermore we have identified the critical values a°_ and a\ at 
which E® = and E\ = respectively. Successively we have calculated the four-, five-, 
and six-body spectra and we have observed the two-level structure. The EFT approach is 
better adapted to describe shallow states and this is confirmed by the close result obtained 
for E\ = -128.8 mK at A = 1 and the LM2M2 result E\ = -127.9 mK, very recently 
published Q. 

The universal character of the structure of these clusters has been studied using Tjon lines, 
that means the relation between E A+1 and E A , and between E A+1 and E\. As illustrated 
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FIG. 10. Relation between E A+1 and E A (upper panel), and between E l A+l and E\ (lower panel) 
for the indicated cases obtained with the TBG-F^biB potential along the (a ,k) plane. 




FIG. 11. (color online). Energies of the states E®, E\, E®, E§, and E® as a function of a -1 
for negative values of the scattering length close to the continuum threshold obtained using the 
TBG potential. The four-particle thresholds are a A _f ~ —19.6 a.u., and aj ~ —39.8 a.u. The 
five-particle thresholds are ~ —12.5 a.u., and a 5 ^ 1 ~ —18.7 a.u. 



in Fig. [TDJ we have obtained an almost linear relation between E\, x and E° A (upper panel) 
and between E\ +1 and E\ (lower panel) in the region from A = 1 to 0.9. As the energy of 
the cluster, E\ or E\, tends to zero the linear relation is lost. 

Since we are describing the lowest bound states, some universal ratios are only approx- 
imately verified, though not very far from the values quoted by other groups in A = 3,4. 
However, in the simple case of TBG we have extended our calculations for A = 4 and A = 5 
up to the four- and five-particle continuum threshold in order to calculate the ratios be- 
tween the thresholds: the values we obtain in the four-body case, (a J /o°)tbg ~ 0.92 and 



{al° /a°_ )tbg ~ 0.45, and in the five-body system (a 5 _l° /a J°) tbg ~ 0.64 and (gl^/gl^tbg. 



0.95, are in accord with the ratios that have been previously predicted 
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5,1 / 4,0\ 



and measured 



27) 



TABLE II. Convergence of the binding energies as a function of the grand angular quantum number 
K using the TBG+H3B potential at A = 0.9, Rq = 10 a.u., and po = 10 a.u., for clusters of A 
particles. We also report the number Ahh of hyperspherical basis elements corresponding to a 
given K. 



K 



iVHH 



E° A (mK) 



E\{mK) 



36 
38 
40 



33649 
42504 
53130 



-166.25945 
-166.25949 
-166.25951 



-6.55041 
-6.79163 
-6.99574 



26 
28 
30 



448800 
724812 
1139544 



-532.75811 
-532.75828 
-532.75834 



-161.96737 
-162.98374 
-163.79689 



18 
20 
22 



709410 
1628328 
3527160 



-1063.8276 
-1063.8311 
-1063.8322 



-513.50956 
-516.42712 
-518.25341 



in literature. 

Another interesting aspect is the uncertainty introduced by the cutoff in the hypercentral 
three-body force. We have observed that with the most natural choice po = R the shallow 
states E\ and E\ result unbound in the last part of the curves. They cross the respective 
threshold E® and E®. Increasing po they result bound again around po ~ 18 a.u.. Increasing 
further p they become again unbound. This last analysis is somehow inconclusive as to 
really understand the cutoff dependence we need to vary both cutoff Rq and po in a coherent 
way; the dependence on the cutoff will eventually reflect the leading order nature of the 
potential we are using, pointing to the necessity of going to a higher order in the EFT 
expansion 19J. Studies along this line are at present under consideration as well as the 



analysis of the two-level spectrum for cluster with more than six particles. 
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